
nc=size(pcmat,1);

hT=25;
nT=(T1-T0)/hT+1;

Tgrid=linspace(T0,T1,nT);

for i=1:nc
    
    figure;
    
    plot(Tgrid,pcmat(i,:),'LineWidth',3);
    hold on;
    
    if joh
        
        plot(Tgrid,pcjmat(i,:,1),'LineWidth',1);
        plot(Tgrid,pcjmat(i,:,2),':','LineWidth',2);
        
    end
    
        
    grid on;
    set(gca,'fontsize',20)
    print('-depsc',['dim' num2str(n) '_c' num2str(i) '.eps']);
    
end